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Abstract. Infrared asymptotic behaviour of a scalar field, passively advected by a 
random shear flow, is studied by means of the field theoretic renormalization group 
and the operator product expansion. The advecting velocity is Gaussian, white in 
time, with correlation function of the form oc 6(t — t')/k d ^ 1+ ^, where k± = \k±\ and 
kj^ is the component of the wave vector, perpendicular to the distinguished direction 
('direction of the flow') — the d-dimensional generalization of the ensemble introduced 
by Avellaneda and Majda [Commun. Math. Phys. 131: 381 (1990)]. The structure 
functions of the scalar field in the infrared range exhibit scaling behaviour with exactly 
known critical dimensions. It is strongly anisotropic in the sense that the dimensions 
related to the directions parallel and perpendicular to the flow are essentially different. 
In contrast to the isotropic Kraichnan's rapid-change model, the structure functions 
show no anomalous (multi) scaling and have finite limits when the integral turbulence 
scale tends to infinity. On the contrary, the dependence of the internal scale (or 
diffusivity coefficient) persists in the infrared range. Generalization to the velocity 
field with a finite correlation time is also obtained. Depending on the relation between 
the exponents in the energy spectrum £ oc k^ 6 and in the dispersion law uj oc 
the infrared behaviour of the model is given by the limits of vanishing or infinite 
correlation time, with the crossover at the ray r\ — 0, e > in the e-r\ plane. The 
physical (Kolmogorov) point e = 8/3, r\ = 4/3 lies inside the domain of stability of the 
rapid-change regime; there is no crossover line going through this point. 
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1. Introduction 

The problem of turbulent advection, being of practical importance in itself, has become 
a cornerstone in studying fully developed hydrodynamical turbulence on the whole pQ. 
On one hand, deviations from the classical Kolmogorov theory — intermittency and 
anomalous scaling [2j |3] — are much stronger pronounced for a passively advected scalar 
field (temperature of the fluid or concentration of impurity) than for the advecting 
turbulent field itself. On the other, the problem of passive advection appears easier 
tractable theoretically. Most remarkable progress was achieved for Kraichnan's rapid- 
change model: for the first time, the anomalous exponents were derived on the basis of 
a dynamical model and within controlled approximations [4, 5j. 

In Kraichnan's model, the turbulent velocity field is modelled by the Gaussian 
distribution with the pair correlation function of the form 

(^•)(xD ^-O^^, (1-1) 

where = 5{j — kikj / k 2 is the transverse projector, k = |k| is the wave number, Dq > 
is an amplitude factor, d is the dimension of the x space and £ is an arbitrary exponent. 
The latter can be viewed as a kind of Holder's exponent, which measures 'roughness' of 
the velocity field; the 'Batchelor limit' £ — > 2 corresponds to smooth velocity, while the 
most realistic (Kolmogorov) value is £ = 4/3 pp. 

The issue of interest is, in particular, the behaviour of the equal-time structure 
functions 

S n (r) = (l9(t,x)-9(t,x')] n ), r=|x-x'| (1.2) 

of the scalar field 6(t,x) in the inertial range £ ^ r ^ C, where I is the dissipation 
length and £ is the integral turbulence scale. Within the so-called zero-mode approach, 
developed in [H [5], it was shown that in the inertial range the functions ( 11.2ft are 
independent of the diffusivity coefficient and have the forms: 

S 2n (r) oc £>-"r"M (r/£) A «, (1.3) 

with negative anomalous exponents A„, whose first terms of the expansions in 1/d pi] 
and £ [5] are the following: 

A n = -2n(n - l)£/d + 0(l/d 2 ) = -2n(n - + 2) + 0(£ 2 ). (1.4) 

Thus the functions ( 11.21) depend on the integral scale and diverge for £ — > oo, in 
contradiction with the classical Kolmogorov theory. 

In [6] and subsequent papers, the field theoretic renormalization group (RG) and 
operator product expansion (OPE) were applied to Kraichnan's model; see £7] for the 
review and references. In the RG approach, the exponent £ plays the part analogous to 
that played by e = 4 — d in Wilson's theory of critical phenomena, while d remains 
a free parameter. The anomalous scaling for the structure functions emerges as a 
consequence of the existence in the corresponding operator product expansions of 
'dangerous' composite fields (composite operators in the field theoretic terminology) 
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of the form (d8) 2n , whose negative critical dimensions are identified with the anomalous 
exponents A n . This allows one to construct a systematic perturbation expansion for 
the anomalous exponents and to calculate them up to the orders £ 2 [6] and £ 3 jH]- 

In this paper, the RG+OPE approach is applied to the model of a passive scalar field 
in a random shear flow: the Gaussian velocity field is oriented along a fixed direction n 
('direction of the flow') and depends only on the coordinates in the subspace orthogonal 
to n. In the momentum space, its correlation function has the form simillar to ( II. lft : 
(vv) oc S(t — t') kj_ d+1 ~^, where k± = |kj_| and kj_ is the component of the momentum k 
perpendicular to n. This model can be viewed as a (^-dimensional generalization of the 
strongly anisotropic velocity ensemble introduced in [9] in connection with the turbulent 
diffusion problem and further studied and generalized in a number of papers [T0]-|19j. 

We show that the inertial-range behaviour of this model appears essentially different 
from the isotropic Kraichnan's model: due to the absence of dangerous composite 
operators, the structure functions (jl.2p have finite limits at C — > oo and thus show no 
anomalous scaling in the sense of ( II. 3ft . On the contrary, dependence on the diffusivity 
(and thus on the dissipation length) persists in the inertial range. Following the 
nomenclature of the monographs [2, 3j, one can say that, in complete contradistinction 
with isotropic Kraichnan's model, the first Kolmogorov hypothesis is valid in the present 
case, while the second hypothesis is violated. 

The paper is organized as follows. The sections EH8] are devoted to the rapid-change 
version of the model (vanishing correlation time); generalization to the finite-correlated 
case is given in section [9j 

In section [2] we give detailed description of the model, present its field theoretic 
formulation and the corresponding diagrammatic technique. In section [3] we analyze 
canonical dimensions and ultraviolet (UV) divergences of the model. We show that, 
after an appropriate extension, the model becomes multiplicatively renormalizable. We 
derive the explicit expression for the only independent renormalization constant, which 
is given exactly by the one-loop approximation. In section H] we derive the differential 
RG equations with exactly known coefficients (/3 function and anomalous dimensions 7) 
and show that they possess an infrared (IR) attractive fixed point, which governs the 
scaling behaviour of the Green functions in the IR range. 

In section [5] we present the corresponding critical dimensions for the basic fields 
and parameters. Our model is strongly anisotropic in the sense that, in contrast to 
previous RG+OPE studies of anisotropic passive advection [20] [22J , it does not include 
parameters that could be tuned to make the velocity statistics isotropic, and hence it 
does not include the isotropic Kraichnan's model as a special case. As an interesting 
consequence, the critical dimensions related to the directions parallel and perpendicular 
to the flow are essentially different. 

Section |6] is devoted to the composite operators. As already mentioned, the key 
role in the RG+OPE approach to anomalous scaling is played by the dimensions of the 
Galilean invariant operators (d6) 2n , built of the scalar gradients [6] [H]- In the isotropic 
case, there is only one such operator for a given n, namely (di8di0) n . In the strongly 
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anisotropic case of a shear flow, there is a set of (n + 1) relevant operators for each n. 
They mix heavily in renormalization and give rise to a set of critical dimensions rather 
than a single A n . Nevertheless, it turns out that exact expressions can be derived for 
these dimensions in our model. Furthermore, in contrast to their counterparts (11 .4p in 
the isotropic Kraichnan's model, they all are positive. 

Sections [7J and [8] apply the results of the preceding analysis to the inertial-range 
asymptotic behaviour of the structure functions (11.21) . In section [7] their behaviour in the 
IR range r 3> £ is established; it turns out that those functions retain the dependence on 
the UV scale £. The inertial range corresponds to the additional condition that r <^ C; 
it is studied by means of the OPE in section [HJ Due to the absence of relevant dangerous 
operators with negative dimensions, the structure functions appear finite for (r/C) — > 
and thus show no anomalous scaling in the sense of ( 11. 31) . The resulting inertial-range 
asymptotic expressions, presenting the main outcome of this study, are summarized in 

(P-jSD- 

In section E3, the generalization of the above results to the velocity ensemble with 
finite correlation time is given. The energy spectrum is taken in the form £ oc k_^ e , 
while the dispersion law is w oc k 2 ^ v . It is shown that the IR behaviour of the model is 
nearly exhausted by the two limiting cases: the rapid-change type behaviour, realized for 
e > i] > (with £ = e — 77 > 0), and the frozen (time- independent) behaviour, realized 
for e > 0, rj < 0. The crossover line between the two regimes is the ray rj = 0, e > in 
the e-T] plane. In contrast to the isotropic case, where the physical (Kolmogorov) point 
e = 8/3, 7] = 4/3 lies exactly on the crossover line between the rapid-change and frozen 
regimes [23j-[26], now this point lies deep inside the domain of stability of the nontrivial 
rapid-change behaviour; there is no crossover line going through this point. This result 
is in agreement with the findings of the exact analysis of the d = (l + l)-dimensional case 
by [THEE] an d in disagreement with [9]-[TT]; this issue is further discussed in section [TOj 
which is also reserved for conclusions. 



2. Description of the model and the field theoretic formulation 

The advection-diffusion equation for the scalar field 9(x) with x = {t,x} has the form 
V t 9 = u d 2 6 + C, (2.1) 

where 

V t = d t + vA (2.2) 

is the Galilean covariant (Lagrangian) derivative, d t = d/dt, di = d/dxi, d 2 = didi is 
the Laplacian, u is the diffusion coefficient and ((t, x) is a Gaussian random noise with 
zero mean and the pair correlation function 

(C(t,x)C(t , ,x'))=5(t-t / )^(r), r = x-x'. (2.3) 

The function C(r) is finite at r = (and we assume the normalization C(0) = 1) and 
rapidly decays for r — > 00; its precise form is inessential. For incompressible fluid, the 
velocity field v = {vi(x)} is transverse due to the continuity relation: diVi = 0. 
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Let n be a unit constant vector that determines some distinguished direction 
('direction of the flow'). Then any vector can be decomposed into the components 
perpendicular and parallel to the flow, for example, x = xj_ + nxn with xi • n = 0. The 
velocity field will be taken in the form 

v = m;(t,x ± ), (2.4) 

where v(t, xj_) is a scalar function independent of xu. Then the incompressibility 
condition is automatically satisfied: 

dm = d\\v(t,x ± ) = 0. (2.5) 

For v(t, x_i_) we assume a Gaussian distribution with zero mean and the pair correlation 
function of the form: 

/dk 
— exp{ik-(x-x')}A,(AO = 

— exp{ik ± ■ (x_L-x! ± )}D v (k ± ), k ± = |kjj (2.6) 

with the scalar coefficient functions 

D v (k) = 27T(5(fc||) D v (k ± ), D v (k ± ) = D kl d+1 -t. (2.7) 

Here and below d is the dimension of the x space, -Do > is a constant amplitude factor 
and £ an arbitrary exponent. The IR regularization in (12 .61) is provided by the cutoff 
k± > m, where m ~ C~ l is the reciprocal of the integral turbulence scale. Its precise 
form is inessential; the sharp cutoff is the most convenient choice from the calculational 
viewpoints. The natural interval for the exponent is < £ < 2, when the so-called 
'effective eddy diffusivity' 

V(rj.) = J T^fcr {1 - exp (ik ± • r ± )} D v (k ± ) (2.8) 

has a finite limit for m — > 0; it includes the most realistic Kolmogorov value £ = 4/3. 

In order to ensure multiplicative renormalizability of the model, it is necessary to 
split the Laplacian in (12. 1 p into the parallel and perpendicular parts d 2 — > dj_ + f d 2 
by introducing a new parameter /o > 0. Here d\ is the Laplacian in the subspace 
orthogonal to the vector n and d\\ = d/dxn. In the anisotropic case, these two terms 
will be renormalized in a different way. Thus equation (12. ip becomes 

V t 6 = u {di + f d 2 }6 + (; (2.9) 

this completes formulation of the model. It remains to note that, for the velocity field 
(12 .4p . the covariant derivative in (12. 2p takes on the form 

V t = d t + v(t, x± )d {l . (2.10) 

Interpretation of the splitting of the Laplacian term in (12. 9 p can be twofold. On one 
hand, stochastic models of the type fl2. ID are phenomenological and, by construction, 
they must include all the IR relevant terms allowed by symmetry. The fact that the 
splitting is required by the renormalization procedure means that it is not forbidden by 
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dimensionality or symmetry considerations and, therefore, it is natural to include the 
general value / 7^ 1 to the model from the very beginning. On the other hand, one 
can insist on studying the original model with /o = 1 and Od covariant Laplacian term, 
although that symmetry is broken to O^-i <8> Z 2 by the interaction with the anisotropic 
velocity ensemble (Z2 is the reflection symmetry x\\ — >■ —x\\). Then the extension of 
the model to the case /o / 1 can be viewed as a purely technical trick which is only 
needed to ensure the multiplicative renormalizability and to derive the RG equations. 
The latter should then be solved with the special initial data corresponding to /o = 1 (in 
renormalized variables this anyway will correspond to general initial data with / ^ 1). 
Since the IR attractive fixed point of the RG equations is unique (see section @|, the 
resulting IR behaviour will be the same as for the general case of the extended model 
with /o^l. 

According to the general theorem (see e.g. chap. 5 of the monograph |27J), our 
stochastic problem is equivalent to the field theoretic model of the extended set of fields 
$ = {9', 9, v} with action functional 

S($) = ^9'D c 9' + 9'{-V t + u (dl + f df { )}e + S v (v), (2.11) 

where D$ is the correlator f]2.3f) . The first few terms represent the De Dominicis- Janssen 
action functional for the stochastic problem (12. ip . (12. 3p at fixed v; it involves auxiliary 
scalar response field 9'{x). All the required integrations over x = {£, x} are implied, 
for example, the coupling term in (12. lip , stemming from the derivative ( I2.10p . in the 
detailed notation has the form: 

-. M ^-/,/, l /MW.(U JVW . (2,2) 

Due to the independence of the velocity field on the longitudinal coordinate Xn, the 
derivative in ( 12. 12ft can also be moved onto the field 9' using integration by parts: 

-9\vd\\)9 = 9(vd\\)9'. (2.13) 

The last term in ( 12. lip corresponds to the Gaussian averaging over v with correlator 
( 12. 6 p and has the form 

S v (v) = ~ J dt J dx ± dx>(i,x x )5- 1 (x ± -x / ± )^(i,x / ± ), (2.14) 

where 

^(rJaD^rf^ (2.15) 

is the kernel of the inverse linear operation D' 1 for the correlation function D v in ( 12. 7p . 

This formulation means that statistical averages of random quantities in the original 
stochastic problem coincide with the Green functions of the field theoretic model with 
action (12. lip , given by functional averages with the weight exp5($). This allows one 
to apply the field theoretic renormalization theory and renormalization group to our 
stochastic problem. 

The action (12. lip corresponds to the Feynman diagrammatic technique with 
three bare propagators: the correlator of the velocity field (vv)o, given by ( 12.61) . 
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I] = 4^ P- + 4^ — he— P Hr- + 






Figure 1. The self-energy operator (I2.20P in the two- loop approximation. 

( 12. 7p . the scalar Green function (in the frequency-momentum and time-momentum 
representations): 

(00') o = {-iu + e(k)}" 1 9(t - t') exp{-e(k)(t - t')} (2.16) 
and the correlator of the scalar field 

c(k) ^ c(k) exp{ _ e(k))t _ t1} _ (217) 



w 2 + e 2 (k) 2e(k) 

Here C(k) is the Fourier transform of the function from (12.31) . e(k) = Va(k\ + /o&|) an d 
6(. . .) is the Heaviside step function, so that the function (12. 16[) is retarded. The only 
vertex (12 . 1 3[) corresponds to the vertex factor 

V(p) = -ip|| = ik\\ , (2.18) 

where p is the momentum argument of the field and k is the momentum of 0'. The 
role of the bare coupling constant (expansion parameter in the ordinary perturbation 
theory) is played by the parameter u> , defined by the relation 

D = w u f , w ~A« (2.19) 

with Dq from (12.71) . The last relation, following from dimensionality considerations, sets 
in the typical UV momentum scale A ~ l/£, the reciprocal of the UV length scale. 

As an example, on figure [1] we show the two-loop approximation of the self-energy 
operator S which enters the Dyson equation 

(0'0)i_ir(w, p) = - {-ioo + vopl + ^o/oPy} + p) (2.20) 

for the 1-irreducible Green function (6'6)i-h- in the frequency-momentum representa- 
tion. The wavy lines denote the velocity correlator (vv)q, while the solid lines correspond 
to the function (96') q; the slashes mark the field 6'. 



3. Canonical dimensions, UV divergences and renormalization 

The analysis of UV divergences is based on the analysis of canonical dimensions 
of the 1-irreducible Green functions. In general, dynamic models have two scales: 
canonical dimension of some quantity F (a field or a parameter in the action functional) 
is completely characterized by two numbers, the frequency dimension dp and the 
momentum dimension dp; see e.g. chap. 5 in [27] . They are determined such that 
[F] ~ [T] ~ d F [L] ~ d p , where L is some length scale and T is the time scale. 
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Our strongly anisotropic model, however, has two independent length scales, related 
to the directions perpendicular and parallel to the vector n, and requires a more detailed 
specification of the canonical dimensions. Namely, one has to introduce two independent 
momentum canonical dimensions dp and d F so that 

[F] ~ [T]-^[Lj-^[L||]- d K 

where L± and L\\ are (independent) length scales in the corresponding subspaces. The 
dimensions are found from the obvious normalization conditions = —d^ ± = 1, 
d\ ± = —di ± = 0, d% ± = d% = 0, d" = —d% = 1, and so on, and from the requirement 
that each term of the action functional (12.111) be dimensionless (with respect to all 
the three independent dimensions separately). The total momentum dimension can be 
found from the relation d F = dp + d F . Then, based on dp and d F , one can introduce 
the total canonical dimension dp = d F + 2d F = d F + d F + 2d F (in the free theory, 
d t oc dj_ oc <9jj), which plays in the theory of renormalization of dynamic models the 
same role as the conventional (momentum) dimension does in static problems. 

The canonical dimensions of the model (12.111) are given in table [TJ including 
renormalized parameters, which will be introduced a bit later. From table [1] it follows 
that our model is logarithmic (the coupling constant Wq ~ is dimensionless) 

at £ = 0, so that the UV divergences manifest themselves as poles in £ in the Green 
functions. 

The total canonical dimension of an arbitrary 1-irreducible Green function T = 
($ . . . $)i_i r is given by the relation 

d r = d + 2 - ^2 = d + 2 - N e ,d e > - N e d e - N v d v . (3.1) 

Here N$> = {Ng, Ngi, N v } are the numbers of corresponding fields entering the function 
T, and the summation over all types of the fields in (13. ip and analogous formulas below 
is always implied. 

Superficial UV divergences, whose removal requires counterterms, can be present 
only in those functions T for which the 'formal index of divergence' dr is a nonnegative 
integer. Dimensional analysis should be augmented by the following observations: 
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(1) In any dynamical model of the type (j2.1ip . 1-irreducible diagrams with Ng/ = 
contain closed contours of retarded propagators (I2.16P and therefore vanish. 

(2) For any 1-irreducible Green function Ng> — Ng = 2Nq, where Nq > is the 
total number of the bare propagators (88) o entering into any of its diagrams. This fact 
is easily checked for any given function; it is illustrated by the function ( 12.201) with 
N e i = Ng = 1 and N = 0; see figure [TJ Obviously, no diagrams with N < can be 
constructed. Therefore, the difference Ng/ — Ng is an even nonnegative integer for any 
nonvanishing function. 

(3) The derivative dn at the vertex 9'vd»9 can be moved onto the field 8' due to 
the transversality of v, see (12.131) . Therefore, in any 1-irreducible diagram it is always 
possible to move the derivative onto any of the external 'tails' 8 or 8', which reduces the 
real index of divergence: d' r = d^ — Ng — Ng>. The fields 8, 8' enter into the counterterms 
only in the form of derivatives d\\8, d\\8'. 

From table [1] and (13. ip we find 

dr = d + 2-(d+l)N e >+N e -N v , <#= (d + 2)(l - ify) - N v . (3.2) 

From these expressions we conclude that for any d, superficial divergences can be present 
only in the 1-irreducible functions (8'8 . . . 6 1 ) i_i r with Ngr = 1 and arbitrary Ng, for which 
dr = 2, d' r = 0. However, all functions with Ng > N$> vanish (see above) and obviously 
do not require counterterms. We are left with the only superficially divergent function 
(0'0)i_i r ; the corresponding counterterm must contain two symbols d\\ and therefore 
reduces to 8'd^ 8. 

Inclusion of this counterterm is reproduced by the multiplicative renormalization 
of the action (12. lip with the only independent renormalization constant Zf. 

= v, fo = fZ f , w = w^Z w , Z w = Zj l . (3.3) 

Here the reference scale \x is an additional parameter of the renormalized theory, u, f 
and w are renormalized analogs of the bare parameters (with the subscript '0') and 
Z = Z(w,^,d) are the renormalization constants. Their relation in (13. 3 p results from 
the absence of renormalization of the contribution with Dq in (12. lip , so that 

D = w iy f = wfihf, (3.4) 

see ( 12 .191) . No renormalization of the fields and the parameter m is required: 

m Q = m, Z m = 1, Z$ = 1 for all $. (3.5) 

Here and below we use the minimal subtraction (MS) scheme, where all renormalization 
constants have the forms '1 + only poles in 

The constant Zf is determined by the requirement that the function (0'6)i-- lt , 
expressed in renormalized variables, be UV finite, that is, finite at £ = 0. We recall that 
the correlator (vv)q contains the 5 function in time, while the propagator (12.161) contains 
the step function. Thus all the multiloop diagrams in the self-energy operator E in (I2.20p 
contain self-contracted chains of the step functions, like e.g. 0(^1—^2)0(^2—^3)0(^3—^1), 
and therefore vanish. (In the frequency representation, all the integrands have the poles 
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in u only in the lower complex half-plane.) This means that the functions £ and (0'6)i_ ir 
are given exactly by the one-loop approximation. 

The analytic expression for the only one-loop diagram has the form 

E(P) = / - / A?, DM ip ^ P ; k \^ (3.6) 
Kyj J 2n J (2rr) d vK ; -vjj + e(p - k) V ; 

with D v from (12. 7p and e(k) from (I2.17p . This expression is independent of the external 
frequency. The prefactor, coming from the vertices (12.181) . can be replaced with — p?, 
due to the presence of the factor S(k^) in (12. 7ft ; this is the diagrammatic analog of the 
relation (I2.13p . Integration over u involves the indeterminacy 

*" 1 8(0), (3.7) 



2ix —iu + e(p — k) 

the step function at the origin; it should be carefully resolved. In our case, the function 
5(t — t') should be understood as the limit of a narrow function which is necessarily 
symmetric in t f-> t', because (12. 6p is a pair correlation function. Thus the quantity in 
(13. 7p must be unambiguously defined by half the sum of the limits: 6(0) = 1/2. 

The integration over kn is trivial due to the presence of the factor 5(k\\) in (12. 7p . 
The result has the form 

w n 2 A) f dk ± 1 

The remaining integration over k^ gives 

E(p) = -pji?om-«^^i (3.9) 

where = 2ir d ^ 2 /T(d/2) with Euler's Gamma function is the surface area of the unit 
sphere in the d- dimensional space. Substituting the expression (13.91) into the Dyson 
equation (12.201) and passing to renormalized parameters with the aid of relations (13. 3 p 
and (13.41) gives 

<M> ,_(«, P) = - u V i + vfZ 1P \} ~P\§^ {^f i (3.10) 

In order to cancel the pole in £, the renormalization constant in the MS scheme has to 
be chosen in the form 

Zf = 1 - j , (3.11) 
where we have absorbed the factor Sd-i/2(27r)( d ~ lS> into the coupling constant. 



4. RG equations and the fixed point 

Let G(e , . . .) be some correlation function in the original model (12. lip and G^e, /1, . . .) 
its analog in the renormalized theory. Here eo is the complete set of bare parameters, e is 
the set of their renormalized counterparts, and the ellipsis stands for the other variables 
like the coordinates/momenta and times/frequences. These functions differ only by 
normalization and the choice of variables, G(eo) = ZcG^e, /i), and can equally be used 
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for the analysis of critical behaviour. (For the correlation functions of the primary fields 
Zq = 1 due to the absence of their renormalization, but it is instructive to discuss a 
more general case.) We use to denote the differential operation /id^ for fixed eo and 
operate on both sides of the last equality with it. This gives the basic RG equation: 

{V RG + lG }G R (e,fi,...) = 0, (4.1) 

where T>rg is the operation D M expressed in the renormalized variables: 

V RG = V„ + f3d w - lf V f . (4.2) 

Here we have written T> x = xd x for any variable x, and the RG functions (the /3 function 
and the anomalous dimensions 7) are defined as 

(3{w)=V fl w = w[-£- lw ], (4.3) 

7e(w) = V^lnZe = (3d u liaZ e for any quantity e. (4.4) 
From the relations (I3.3p -( l33|) it follows 

lv = 7m = 7* = and 7/ = -j w , (4.5) 

while from (14. 31) and ( I4.4j) for e = w one easily derives 

-£V w lnZ w -£w 

1 + V w in Z w 1 + V w In Z w 

Substituting (13.111) into (14. 6p one obtains exact expressions 

7/ = -Jw = w, P = w[-£ + w]. (4.7) 

It is well known that IR asymptotic behaviour of the Green functions is governed by 
IR attractive fixed points of the RG equations, defined by the relations = and 

> 0. From (14.71) it follows that our model has a fixed point 

w* = e, P\w*) = £ (4.8) 
which is positive and IR attractive for ^ > 0. At this point 

1} = = £• (4-9) 
Here and below we denote 7* = j e (w*). We also stress that all the expressions in (14.81) 
and (14. 9p are exact. 

5. Critical scaling and critical dimensions 

In the leading order of the IR asymptotic behaviour the Green functions satisfy the RG 
equation with the substitution w — > w*, which gives 

{X? m - 7 P/ + 7g} G R (e,/i,...) = 0. (5.1) 

Canonical scale invariance of the function G R with respect to the three independent 
canonical dimensions (see section [3]) can be expressed by the differential equations of 
the form 

d* F V F -d* G I G fl (e,/i,...) =0, (5.2) 
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where the sum runs over all arguments of the Green function, including the 
coordinates/momenta and times/frequences, and d* F = dp, dp and d F are the 
corresponding canonical dimensions. In the time-coordinate representation 

d* G = Yl N * d * = Ne ' d v + N ° d *o + N ^ ( 5 - 3 ) 

where Nq> are the numbers of the fields entering the Green function, cf. equation (13. ip . 
From table [1] we find 

{-V t + V v -d" G } G R (e,/i,...) =0, 

{-V x +V fi -2V„ + 2V f + V m -d G } G R (e, »,...) = 0, 

-2V f - 4} G R (e, //,...) =0, (5.4) 

where for definiteness we use the time-coordinate representation and denote V± = 
x±d/dx±, V\\ = x\\d/dx\\. 

The equations of the type (15. ip . ( 15. 2ft and ( 15. 4p describe the scaling behaviour of the 
function G R upon the dilation of a part of its parameters: a parameter is dilated if the 
corresponding derivative enters the equation; otherwise it is kept fixed. We are interested 
in the IR scaling behaviour, in which all the IR relevant parameters (coordinates, 
times, integral scale) are dilated, while the irrelevant parameters (diffusivity coefficients, 
coupling constant) are fixed. Thus we combine the equations ( 15. ip and ( 15. 4 p so that the 
derivatives with respect to the IR irrelevant parameters be eliminated; this gives the 
desired equation which describes the IR scaling behaviour: 

{V ± + A||P|| + A W V U + A m V m - A G } G R = 0. (5.5) 

Here = 1 is the normalization condition, while the critical dimensions of any other 
IR relevant parameter F is given by the general expression 

A F = 4 + A||4 + A aJ ^ + 7 ^, (5.6) 

where 

A w = 2- 7 :, A,| = (2 + 7 ;)/2. (5.7) 

Then using ( 14. 5p . ( 14. 9 p and the data from table [1] we obtain the following exact 
expressions for the dimensions: 

A w = 2, A,, = 1 + £/2, A ei = d + 1 + e/2, A = -1, A m = 1. (5.8) 

Solution of the RG equations for the structure functions ( 11 ,2p will be discussed in 
section 

6. Critical dimensions of the composite operators 

The key role in the following will be played by the critical dimensions of certain 
composite fields ("composite operators" in quantum-field terminology). Detailed 
exposition of the renormalization procedure of composite operators can be found in 
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[27]. In general, counterterms to a given operator F are determined by all possible 1- 
irreducible Green functions with one operator F and arbitrary number of primary fields, 
r = (F(x)^(xi) . . . $(x2))i-ir- The total canonical dimension of this function (formal 
index of divergence) is 



where dp is the dimension of the operator and the summation over all types of fields is 
implied, cf. expression (13. ip . For superficially divergent diagrams dr is a nonnegative 
integer. 

We begin with the simplest operators 8 n (x), which enter the structure functions 
(11. 2p . Using the data in tabled] we obtain dp = ndg = —n and dr = —n + Ng — N v — 
(d + l)Ng>. From the analysis of the diagrams it follows that the total number of the 
fields 9 entering the function V cannot exceed the number of the fields 9 in the operator 
9 n itself, that is, Ng < n. This is a direct consequence of the linearity of the original 
stochastic equation ( 12. ip in the field 9: the solid lines in the diagrams cannot branch. 
Therefore, the superficial divergence can only exist in the functions with N v = Ng> = 
and arbitrary value of n = Ng, for which the formal index vanishes, dr = 0. However, in 
any nontrivial diagram at least one of Ng external 'tails' of the field 9 is attached to a 
vertex 9'(vd)9, at least one derivative d appears as an extra factor in the diagram, and, 
consequently, the real index of divergence d' T is necessarily negative. 

This means that the operator 9 n is in fact UV finite and requires no counterterms 
at all: 9 n = Z [9 tl ]r with Z — 1. It then follows that its critical dimension, given by 
the expression (I5.6P without the contribution of 7^, is simply given by the sum of the 
critical dimensions of its constituents: 



In Kraichnan's rapid-change model, the anomalous exponents (11.41) are identified 
with the critical dimensions of the scalar composite operators of the form (d9) 2n , the 
operators with minimal canonical dimension (namely, dp = 0) that are invariant with 
respect to the shift 9 —> 9 + const; see [6] [7]. For that isotropic case, the scalar operator 
of the needed form is unique for any given n: F n = (di9di9) n . The operator F n does not 
admix in renormalization to F^ with k < n, the corresponding renormalization matrix 
is triangular, and the dimensions A„ are determined by its diagonal elements. (The 
operators F n can be treated as if they were multiplicatively renormalizable.) For the 
existence of the singular behaviour of the structure functions ( 11.31) at (r /£)—>• it is 
crucial that the dimensions A n are negative. 

In our case the Od isotropy is broken to Od-i <S> Z2, and for any n one can construct 
a set of (n + 1) different operators of the form (d9) 2n , invariant under the residual 
symmetry, namely: 




(6.1) 



A[9 n ] 



nAg = —n. 



(6.2) 



F K 



s 



(di9di9) k (d\\9d\\9)\ k + s = n. 



(6.3) 
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r = • + 




Figure 2. The one- loop and a part of the two-loop contributions to the generating 
functional (16.61). 



Here 8^ is the derivative in the subspace, d\\ = d/dx\\ and the summation over the 
vector index i — 1, . . . , d — 1 is implied. In particular, for n — 1 such a set includes two 
operators 

F 1)0 = (dtBdtO), F 0A = (d l{ 9) 2 , 
— three operators 



(6.4) 



for n 



2,0 



(6.5) 



= and 
following 



{Dudley, F hl = {diedte){d\\ef, f 0j2 = (d^, 

and so on. 

For all operators (I6.3p . from table [1] and equation (16. ip we find dp 
dr = Ng — N v — (d + l)Ng>, with the necessary condition Ng <2n = 2(k + s) 
from the linearity of the model (cf. the discussion for 9 n above). From the form of the 
vertex 9'(vd)9 and the operators Fk )S themselves it follows that the fields 9, 9' enter the 
counterterms only in the form of derivatives, 89, 89', so that the real index of divergence 
is d' T = d r — N e — Ng> = —N v — (d + 2)Ng>. Thus the superficial divergences can only 
exist in the Green functions with N v = N& = and arbitrary Ng < 2n (then d' r = 0), 
and the corresponding operator counterterm necessarily reduces to the form F m ^i with 
(m + l) < (k + s). 

Therefore, the operators of the type Ff. a can mix only with each other in 
renormalization, the corresponding infinite renormalization matrix Z = {Zk tS - mt i} is 
in fact block-triangular, Zk jS ;m,l — for (m + l) > (k + s), and the critical dimensions 
associated with the set of operators with given n = (k + s) are determined by its finite 
diagonal blocks Z n = {Zk, S ;m,i} with the fixed value of the sum n = (m + I) = (k + s). 

In the following, we will show that the counterterms and the renormalization 
matrices for these operators are fully specified by the one-loop approximations. 

It is convenient to work with the generating functional of 1-irreducible correlation 
functions with one operator F^ s and arbitrary number of the primary fields 9: 



k,s 



( x - 9) = f dxi ■ ■ ■ f dx p (F Ks (x)9(x 1 ) ■ ■ ■ 9(x p )\ _ 9(x x ) ■ ■ ■ 9(x p ). 

(6.6) 

In the diagrammatic expansion of the functional r^, the leading (loopless) 
approximation is given by the operator Fk s (x) itself. The only one-loop diagram and 



Inertial-range behaviour of a passive scalar field in a random shear flow 



15 



two of the five two-loop diagrams are shown in figure [21 The diagrammatic notation 
was explained in section [3] in connection with figure [TJ the new elements are the thick 
dots that denote the vertex factors of the composite field Fky. 

vP(x- xi x,) - 5lFk >°W ( 6 7 ) 

v kA x > x i>---> x v se(x 1 )---se(x l y [ } 

where I is the number of (88')o-\mes attached to the vertex (16. 7ft . The variational 
derivative acts on such operators as follows: 
5 



-F ktS (x) = 2kF k . hs (x) df9(x) df5(x - x') + 



56{x') 

+ 2s F fc , s _i(x) d\\6(x)d\\5(x - x'), (6.8) 

so that the vertex factor (16.71) is in fact local. 

The derivatives dn at the external vertices (two lower vertices in the first two 
diagrams on figure [2] and three lower vertices in the last diagram) correspond to external 
momenta; they act onto the fields 9 attached to those vertices and form the external 
factors ~ (d\0. The remaining expressions diverge logarithmically and therefore can 
be calculated at all external momenta equal to zero; the IR regularization is provided 
by the cutoff m in the velocity correlator (12.71) . 

It is easy to see that the number of independent loops in any diagram of the 
functional (16.61) equals to the number of wavy lines; thus the independent integration 
momenta (denoted as k and q below) can always be assigned to the velocity correlators. 
Then the integrals corresponding to the diagrams in figure [2] take on the forms: 

Jk f<L ' hh ' J [D v {k){d\\0{x)f] (6.9) 



(2vr) d J 2vr(a; 2 + e 2 (k)) 
for the one-loop diagram, 



(2vr) d J (27r) d J 2tt J 2tt 

(fc + g)t(fc + g)j 

X {oo 2 + e 2 (k + q) + e 2 (q)} {Q 2 + e 2 (q)} 

for the first two-loop diagram and 

rfk f rfq f duj f dfl 



(6.10) 



(2ir) d J (2ir) d J 2tt J 2tt 

ikiqj(k + q)i 



[D v {k){i % d\\9{x))\ [D v {q){d\\9{x)f] x 

(6.11) 



{Q? + e 2 (q)} {-iu + e(k)} {i(w + SI) + e(k + q)} 

for the last one. Here e(k) = v(k\ + fk?<) is the renormalized analog of the quantity 
from (ETTjl . 

The momenta with free vector indices stem from the vertex factor (16. 7p and 
correspond to the derivatives acting onto the spatial 5 functions in expressions like 
( 16. 8ft . Each factor in square brackets comes from one velocity correlator (12. 6p and two 
vertices (I2.18P attached to it; D v being the function (12. 7p . Thus the integrand for any 
diagram contains the factor D v (k) ~ 5(k\\) for each independent integration momentum 
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k. On the other hand, any diagram of the functional (16. 6p with more than one loop 
necessarily contains at least one internal vertex (j2.18p . the corresponding vertex factor 
is proportional to certain linear combination of independent parallel momenta k\\, q\\ and 
so on, and the integrand as a whole vanishes due to the presence of the corresponding 
5 functions (we recall that all the external momenta are set equal to zero). Thus the 
only contribution to the counterterm comes from the one-loop diagram. For the same 
reason, the integral (16. 9p is nontrivial only if the both momenta in the factor kikj are 
perpendicular to the vector n; otherwise ki = nikn and the above mechanism works. So 
one can replace fcj — > kj- in (16.91) . and in the vertex keep only the term 

Thus the operator F^s in the divergent part of the one-loop diagram loses two factors 
d x 8 in (I6.12p and acquires two factors dnd in (I6.9p . so that the counterterm to the 
functional (16. 6p reduces to the form We conclude that the renormalized 

operator corresponding to Fk tS has the form F^ s = Fk tS + aFk-i tS +i- Resolving these 
relations for F^ s gives 

l=k 

^M = ^M + E a ^M^ ( 6 - 13 ) 
1=1 

with some coefficients ai and a\ = —a. Thus the renormalization matrix Z n = {Zk tS - m j} 
in the matrix relation F = Z n F R is triangular (with a natural numbering of operators 
in the family with fixed k + s = n) and its diagonal elements are all equal to unity, 
Zk, s -,k,s — 1- The other nonvanishing elements Zk )S -,k-i,s+i are determined by the one- 
loop diagram (I6.9p . but we will not need their explicit expressions. Indeed, the matrix 
of critical dimensions for the operators (16.31) is given by the expressions (I5.6p - (I5.8I) . in 
which d* F should be understood as diagonal matrices of canonical dimensions of the 
operators Fk, s and 7^ = Z^V^Zn is the matrix of their anomalous dimensions. The 
latter is triangular with vanishing diagonal elements 7jt jS; jt jS = and nontrivial , jk,s-,k-i,s+i 
with k — 1, . . . , I. Thus the critical dimensions Afc jS related to the set Fk, s , given by the 
eigenvalues of the matrix (15. 6p . coincide with its diagonal elements. From table [TJ for 
the operators we find 

dp = 2k, d ll F = 2s, dp = -(k + s), 

so that 

A M = 2k + 2sA|| — (k + s)A ul 

with no contribution from 7^, which along with expressions (15. 8 p gives the final exact 
result 

A M = s£ (6.14) 

for the critical dimensions related to the operators (16. 3p . In contrast to the result (II. 4ft 
for the isotropic Kraichnan's model, they have no corrections of order 0(£ 2 ) and higher 
and are positive for all k, s and £ > 0. 
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7. RG and the IR scaling behaviour of the structure functions 

The results of the preceding sections allow one to find the IR scaling behaviour of the 
correlation functions in our model. For generality, consider the different-time structure 
functions 

5 2 „(r,r ± ,r||) = ([^,x)-^',x')] 2n ), (7.1) 

where t = t' — t, r± — |xj_ — x'_J and r\\ = \x\\ — x',,\. The equal-time functions (11. 2p are 
obtained for r = 0. 

The function (17.11) is a linear combination of the two-point correlators 
(6 k (t, x.)6 s (t', x')) with the fixed k + s = 2n. Due to simple exact relations (16.21) for 
the dimensions of the operators 9 n , the critical dimensions of these correlators are all 
equal, A k + A s = —(k + s) = — 2n = 2nA 9 , and the function S 2 n in the IR range behaves 
as a single object. Its IR scaling form can be easily written using the IR scaling equation 
(I5.5P with known critical dimensions (15.81) of the basic IR relevant parameters, but it is 
instructive to discuss its derivation from the RG equation in more detail. 

From the dimensional considerations one can write 

S 2n = v~ n r^T (fir ± , w, r\\/f l ' 2 r^ urr 2 ± , mr x ) , (7.2) 

where J-"(. . .) are some functions of completely dimensionless arguments (that is, 
dimensionless with respect to all the canonical dimensions d* F separately; see section [3]). 

Since the operators 9 n are not renormalized, the function S 2n satisfies the 
homogeneous RG equation ( 14.1 ft with 7g = 0. Its solution can be represented in the 
form 

S 2n = u^rfF (l,w, n /f l / 2 r x ,I?T/ri,mr ± ) , (7.3) 

with the same functions J 7 . The invariant variables e are solutions of the homogeneous 
equation (14. ip normalized with respect to e at fir± = 1. They can be expressed in the 
original (bare) parameters eo using the relations 

w = wr~^Z w {w), f = fZ f (w), v = v , rh = m , (7.4) 

with the renormalization constants from (I3.3P and (13. 5p . The representations (17. 3p and 
(I7.4p are valid because both their sides satisfy the same RG equation and coincide in the 
normalization point fir± = 1 due to the definition of e and the relations (13. 3p . The simple 
result for v and rh follows from the fact that these parameters are not renormalized. 

The advantage of the representation (17. 3p is that the invariant variables have simple 
asymptotic behaviour for Ar± — > oo (or fir± — > oo in renormalized variables). The exact 
explicit expression for the invariant charge 

(7.5) 

is derived from the relations (13.31) . (13.111) and (I7.4p . From (17. 5p it follows that, for £ > 
and rj_ — > oo, the invariant charge tends to the IR attractive fixed point of the RG 



w 



1 + w^rj/wo 
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equation: w — > W* = £. For /, eliminating the renormalization constant Z w = Zj 1 from 
(J73D gives 

/ = fow r{/w -> f w r{/w*. (7.6) 

Substituting these expressions into (17.31) gives the desired asymptotic expression for the 
structure functions in the IR range Ar_|_ 1: 

S 2n = v^ n r^7 (l, w*, r||/(/ ^o) 1/2 ri +5/2 , u r/rl, mr ± ) . (7.7) 

It can be made more transparent by discarding in the notation all the IR irrelevant 
parameters: 

S?, n = r^TZ ^r\\r~ 1 ~ i/2 ,r/r 2 ± ,mr ± ^ = 

= rl 2nAe H U\/r^\r/r^,mr^A (7.8) 

with certain scaling functions TZ. Of course, the last expression, in which we have 
used the explicit forms (15.81) of the critical dimensions, could be derived directly 
from the general IR scaling equation (I5.5p . but the more detailed representation (17.71) 
gives additional interesting information about the dependence on the IR irrelevant 
papameters. 

In the isotropic Kraichnan's model, the structure functions in the IR range Ar± ^> 1 
depend only on the IR (integral) scale £ = 1/m and the amplitude D Q = wo^o entering 
the velocity correlator (II. ip . but not on the diffusivity coefficient z/ and the coupling 
constant Wo ~ separately; see (11.31) . The same effect takes place for the stochastic 
Navier-Stokes equation; see e.g. the discussion in chapter 6 of [27J. This fact is 
in agreement with the second Kolmogorov hypothesis about the independence of the 
correlation functions in the IR range of the parameters, related to the UV (dissipation) 
scale: viscosity coefficient for the velocity and diffusivity coefficient for the passive scalar 
field; see e.g. [2]. In the case at hand, the UV parameters z/ , fo and wo ~ survive 
in the IR asymptotic expression (17. 7p for the structure functions (they do not form the 
combination D = wo^o/o even if we set fo = 1). Thus we may conclude that, in contrast 
to the isotropic case, the second Kolmogorov hypothesis is invalid for the shear flow. 

8. OPE and the inertial-range behaviour of the structure functions 

The asymptotic representations (17.71) and (17.81) hold in the IR asymtotic range, specified 
by the inequality Ar± ^> 1 (or fir± 3> 1 in renormalized variables), while the other 
arguments of the scaling functions 1Z are kept finite. Inertial range corresponds to 
the additional condition mr^ 1. The form of the scaling functions J 7 or TZ is not 
determined by the RG equations alone. In order to study the limit mr± — > in the 
structure functions, one should combine the plain RG with the OPE techniques; see 
[6]-[7] for Kraichnan's model. For simplicity, below we concentrate on the equal-time 
functions, setting r = in (17.71) and (17. 8p . 
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According to the OPE, the behaviour of the quantities entering into the structure 
functions (11.2j) for r = x — x' — > and fixed x + x' is given by the infinite sum 

[6>(t,x) -6(t,*T = ]TCMr)F ft, ^) , (8.1) 

F 

where Cf are coefficients regular in m 2 and F are all possible renormalized local 
composite operators allowed by the symmetry. More precisely, the operators entering 
the OPE are those which appear in the plain Taylor expansion, and all the operators 
that admix to them in renormalization. 

The structure functions (11.21) are obtained by averaging equation (18.1 j) with the 
weight exp Sr, where Sr is the renormalized action; the mean values (F) appear on the 
right-hand side. Their asymptotic behaviour for m — > is found from the corresponding 
RG equations and has the form (F) oc m Ap . Then substituting the OPE into the 
asymptotic expression (17. 8p gives 

S2n(r) = rl n J2^F (rii/r^mrj.) (mr ± ) A ^, (8.2) 

F 

with coefficients Ap regular in (mr±) 2 . 

Due to the linearity of the model, the number of the fields 9 in such operators 
F cannot exceed their number in the left-hand side. Due to the invariance of the 
action functional (12. lip and the quantity in (18. ip with respect to the shift of the field, 
9 — > 9 + const, the operators entering the OPE must also obey this symmetry. It 
can always be assumed that the expansion (18.11) is made in irreducible tensors (scalars, 
vectors and traceless tensors); owing to the symmetries of our model, only contributions 
of the scalar operators survive in (18. 2p . The leading contributions are determined by 
the operators with minimal critical dimensions Ap = dp + 0(£). 

It then follows that the leading terms of the small-mr^ behaviour of the expression 
(18 .2p for the function Sin are given by the scalar operators Fk yS from (I6.3P with k + s < n. 
Their dimensions Ak yS = in (16.141) are all nonnegative; the leading term is given by 
the operators with s = (including the simplest F — 1), the operators with s > 1 
determine the corrections vanishing for mr± — > 0. 

We conclude that the function S^n remains finite at m — 0. Thus in the inertial 
range and for r = expression (17. 8p becomes 

S 2n = r- ± 2nA °n(r {l /r^), (8.3) 
which for r» = turns to the simple power law 

S 2n ~ r~ 2n ^ = rt, (8.4) 
while for r± = one obtains 

„ -2nAfl/A|| 2n/(l+£/2) / Q r\ 

t>2n ~ r,, =r \\ \ ' 

Thus the inertial-range behaviour in the rapid-change model of a shear flow differs 
essentially from that in the isotropic model. In contrast to the behaviour (11.31) . (11.41) 
of the latter, singular at mr — > 0, there is no anomalous (multi) scaling in the case at 
hand. 
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9. Finite correlation time 

In this section we will consider the case of the Gaussian velocity field with a finite 
correlation time. Now the function D v in the correlator (12.61) . (12. 7p depends on frequency 
u and will be chosen in the form 

DJu.k]) — — t, . 9.1) 

u 2 + [u is k 1 _ '\ z 

The function (19. ip involves two independent exponents r] and e, which in the RG 
approach play the role of two formal expansion parameters. The former defines the 
dispersion law oo(k±) ~ k 2 _^ v (z = 2 — r] in the notation of [9]-[T9]), while the latter 
governs the behaviour of the one- dimensional velocity spectrum, related to the equal- 
time correlator: 

£(k ± ) ~ k d ± ~ 2 [ pD v (u, k ± ) = (g f ^/2u )k]_- £ ; (9.2) 

./ 27T 



this explains the choice of the exponent in the numerator of ( 19. ip . 

The correlator (19. ip involves two important special cases, which, as we will see, 
nearly exhaust possible IR behaviour of the model. The limit uq — > at fixed g' Q = go/uo 
corresponds to the case of time-independent ('frozen') velocity field, when (19. ip turns to 
D v ~ 5{u))k^ d ~ e . The limit u — >■ oo at fixed = go/u^ returns us to the rapid-change 
case (12. 6p . (12 .7p with £ = e — rj and w = g^. 

The role of coupling constants will be played by the parameters (see below) 

So-A^-A^, u ~A\ ^~A £ , <tf~A*, (9.3) 

where A is a typical UV momentum scale, cf. (I2.19p . The model is logarithmic (the 
couplings in (19.31) are all dimensionless) at r\ = e = 0, so that the UV divergences 
manifest themselves as poles in t], e and their linear combinations. 

As a rule, synthetic velocity ensembles with a finite correlation time suffer from 
the lack of Galilean invariance, which can lead to some physical pathologies; see e.g. 
the discussion in [28]. Surprisingly enough, in our strongly anisotropic case the action 
functional (12.111) with the correlator ( 19. II) in (12.141) is invariant with respect to the 
Galilean transformation of the fields 

0(t, x) ->- 0(i, x + ut), 6'(t, x) ->- 6'(t, x + ut), 

v(t, x) -)• v(t, x + ut) - u, (9.4) 

where the transformation parameter has the form u = nu with vector n from (12.41) . so 
that the scalar coefficient in ( 12.41) changes as v (t, xj_) — > f(t,xj_) — u and the arguments 
xi of all the fields in (19.41) remain intact. This fact can be interpreted as follows. 
Consider the stochastic Navier-Stokes equation 

d t Vi + (vidi)vi + dip = R(d)vi + <p i} (9.5) 

where R(d) is some linear differential operation, p = —d~ 2 (diVi)(diVi) is the pressure and 
(f>i is a white-in-time random force. The equation (19.51) is of course Galilean covariant. 
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For the velocity field of the form ( 12. 4 p nonlinear terms in (19. 5p vanish due to the 
independence of the scalar coefficient v on sen: v^d^Vi = Uivd\\v = and similarly for 
the pressure. Thus the equation ( 19. 5 p becomes in fact linear and generates a Gaussian 
velocity field. Its correlator coincides with (19. ip if one choses R(k) = u^v^k 2 ^ and 
0, = <M with (#) = g ulf 5{t - f ) 5{k\\)k"'- d ~ {e+v) . 

The analysis similar to that performed in section [3] shows that the UV divergences in 
the model are removed by the only counterterm 9'd 2 9. Thus the model is multiplicatively 
renormalizable with the only independent renormalization constant Zf. 

9o = g^ £+r 'Zg, u = u/j, v Z u , f = fZ f , u = vZ„, (9.6) 

where 

Z g = Zj\ Z U = Z U = 1. (9.7) 

The constant Zf is found from the requirement that the 1-irreducible Green function 
(#'#)i_i r , expressed in renormalized variables, be finite at e = 77 = 0. Since the 
counterterm has the form 9'd 2 9, the UV divergent part of the self-energy operator 
£(u;, p) in (I2.20p is proportional to p| and it is sufficient to calculate it at vanishing 
external frequency u — 0. Furthermore, pjj is isolated in any of its diagram as an 
extra factor due to the two external vertices (12.181) : see figure [1] and discussion below 
equation ( 13.61) in section [31 Thus in the rest of the corresponding integrals one can set 
p = 0, and the mechanism described below equations (I6.9I) - (I6.11I) 'kills' all the diagrams 
with more than one loop. We stress that, in contrast to the rapid-change case, these 
diagrams are nontrivial; it is only their divergent parts that vanish. 

We are left with the only one-loop diagram; the corresponding frequency integral 
is well-defined, and we obtain 

fduf d^ g^^fk 5 ± ' d -^ 1 

[P) P " J 27T J k±>m (27T)*" 1 ^ + {U^vk 2 -y -ICO + Vk\ - 

- P H ~^lv~ J k±>n (2^=1 1 + «W ' (,) 

The integral in ( 19. 8 p and thus the renormalization constant Zf can be calculated as 

expansions in u, the individual terms would contain the poles ~ l/(e + sr)) with s > 1; 

cf. equations (3.16), (3.17) in [21] and (3.18)-(3.20) in [25]. The calculation can be 

simplified [26] by observing that the anomalous dimensions in the MS scheme are 

independent of the UV regulators like e and 77 (at least in the one-loop approximation for 

our model with several such regulators), so that it is sufficient to calculate the constant 

Zf for 77 = 0, when the integral ( 19. 8ft remains finite. This gives 

2 gifv 2 f f du f dk A 
E(p)=p,| 



and 



2u(u+l)J 2tt J k±>m (2Tr) d -i k d - l+e 

2nt.ii -hi) (2-)''- J V/;;, 1 



P i 9 ; 2f ^ Sd ~\ (^y 1 - (9.9) 



(a 1 c\\ 



Inertial-range behaviour of a passive scalar field in a random shear flow 



22 



where the factor Sd_i/2(27r) <i ~ 1 is absorbed into g. We stress that the expression for 
7/ is exact: in contrast to the isotropic case [26J it has no corrections of order g 2 and 
higher. 

In general, coordinates g^ of the fixed points in a problem with several coupling 
constants are found from the requirement that the /3-functions, corresponding to all 
renormalized couplings vanish. The type of a fixed point is determined by the matrix 
Q with the elements = dPi/dg^, where Pi is the full set of /3-functions and g^ is the 
full set of couplings. For an IR attractive fixed point the matrix Q is positive, that is, 
the real parts of all its eigenvalues are positive. In our case gi = {g, u}\ although u is not 
an expansion parameter in the perturbation theory, the renormalization constants and 
anomalous dimensions depend on it, and it should be treated as an additional coupling 
constant. The corresponding functions (3 (gi) = Tt^gi have the forms: 

Pg = 9[-{ £ + v) +7/], P U = -UT), j f = g/u(u+l). (9.11) 

Since d/3 u /dg = 0, the matrix Q is triangular and its eigenvalues are simply given by 
the diagonal elements fl g and Q u . 

Analysis of the P functions (19. lip reveals several possible fixed points: 

1) «* = with Q u = —rj. Obviously, this corresponds to the frozen case (see the 
remark below equation (19. 2p . and it is convenient to pass from g to the new coupling 
g' = g/u with the (3 function P g > = V^g' = {l/u)P g — (g/u 2 )(3 u , which remains finite at 
u = 0: j5 g i = g'[—e + 7/] = g'\—e + g'\. Thus we find two fixed points: 

1°) 9* — 0; u * — with Q g = —e, Q u = —r] and 7^ = 0, IR attractive for e < 0, 
r] < 0; 

1^) 9* = £ 5 u * = with Qg = 7j = e, Q u = —v, IR attractive for e > 0, i] < 0. 

2) u* = 00. This corresponds to the rapid-change case (12. 6p . (12. 7p . and it 
is convenient to pass to new charges v = 1/u, g" = g/u 2 with the /3 functions 
P v = -(l/u 2 )P u = vi] and Pg,, = (l/u 2 )P g - (2g/u s )P u = + 7/ ] with £ = e-ri, 
which for v = gives fi g n = g"[— £ + g"\. 

Thus we find two more fixed points: 

2a) g'l = 0, t>* = with Q g = fl u = r] and 7^ = 0, IR attractive for £ < (that 
is, e < 77), 77 > 0; 

26) 9* — e i v * — with fig = 7^ = £, Q u = r], IR attractive for £ > (that is, 

6 > Tj), T] > 0. 

For the special case 77 = the function /3 U and the eigenvalue fl u vanish identically, 
and the nontrivial fixed point, attractive for e > 0, becomes degenerate: ^*/u*(u* + l) = 
e. However, the anomalous dimension 7^ = e is independent of its coordinate. 

In figure [3] we show the domains in the E-rj plane, where the fixed points listed 
above are IR attractive. The boundaries of all domains are given by straight lines; there 
are neither gaps nor overlaps between the domains. This fact is exact due to the absence 




Figure 3. Domains of IR stability of the fixed points in the model (19.11) . The numbers 
in boxes correspond to the fixed points (la)— (26) in the text. P is the Kolmogorov 
point e = 2rj — 8/3. 



of the higher-order terms in the functions (19. lip . We also stress that the Kolmogorov 
values of the exponents e = 8/3, r\ = 4/3 lie deep inside the domain of stability of the 
nontrivial rapid-change point (26); there is no borderline going through this point. 

The RG and OPE analysis of the preceding sections equally applies to the model 
(19. ip with a finite correlation time. The IR asymptotic expressions for the structure 
functions in the nontrivial regime (26) have the forms ( I8.2p -( l8~5l) with dimensions (I5.8p . 
(16. 2 p and (I6.14p . For the regime (16) one has to replace £ with e: 

A w = 2, A|| = l + e/2, A e , = d + l + e/2, A m = 1, 

A[0 n ] = nA e = -n, A M = se. (9.12) 

For the trivial regimes (a), the dimensions coincide with their canonical values: A w = 2, 
A || = 1 and so on. In these regimes, turbulent advection is irrelevant in the leading order 
of the IR behaviour, and the difference between the frozen (la) and the rapid-change (2a) 
cases manifests itself only in correction terms (for example, in UV corrections governed 
by the eigenvalues of the matrix Q). Probably for this reason the difference between such 
regimes was not mentioned, e.g., in [9] [TO], [IE], [E]. Adopting the terminology of phase 
transitions, used in those papers, we can say that a first-order transition occurs when 
the point in the e-r\ plane, representing the state of the system, moves continuously 
from domain (16) to (26) and crosses the line rj = 0, e > 0: the effective correlation time 
jumps discontinuously from infinity to zero. On the other hand, if the correlation time 
changes continuously and passes all finite values from infinity to zero, the point in the 
e—r\ plane 'stacks' on the line rj = 0, e > (we recall that the value of is arbitrary 
for the special fixed point with 77 = 0). It remains to note that the critical dimensions 
change continuously from ( 19 . 1 2[) to the 'rapid-change' values (15. 8p . (16. 2p and (I6.14p . 
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10. Conclusion 

We studied inertial-range behaviour of a passive scalar in a random shear flow, modelled 
by a (i-dimensional generalization of the Gaussian ensemble introduced in [9]. The case 
of vanishing correlation time ( 12.4l) -( |2~7l) is discussed in detail, but all the results are 
generalized to the case of finite correlation time (19.1 ft with the spectrum E oc k]^ £ and 
the dispersion law uj ~ k^ v . 

It turns out that possible nontrivial types of the IR behaviour reduce to the two 
limiting cases: the rapid-change type behaviour, realized e > rj > 0, and the frozen 
(time-independent) behaviour, realized for e > 0, 77 < 0. The structure functions in the 
IR range exhibit scaling behaviour of the form (17. 8p . and the corresponding dimensions 
are found exactly: (15. 8 p for the rapid-change case and (I9.12p for the frozen case. 

The resulting inertial-range asymptotic expressions, presenting the main outcome 
of this study, are summarized in (18 . 3 [) — ( T8~5|) . 

In a few respects, the IR behaviour of the model differs drastically from that of the 
isotropic Kraichnan's rapid-change model: 

(i) The scaling is strongly anisotropic in the sense that the critical dimensions, 
related to the directions parallel and perpendicular to the flow, are different. 

(ii) The structure functions in the IR range r 3> i retain the dependence on the 
UV scale t (the second Kolmogorov hypothesis is violated). 

(iii) Due to the absence of relevant dangerous operators with negative dimensions 
in the corresponding OPE, the structure functions appear finite for r ^ C, where £ 
is the integral scale, and thus reveal no anomalous ( mult i) scaling in the sense of (II. 3p . 
Thus the first Kolmogorov hypothesis is valid for our model. 

(iv) For the finite-correlated isotropic case, the Kolmogorov values e/2 = rj = 4/3 
lie exactly on the crossover line between the rapid-change and frozen regimes |23]-|26|. 
For the present model, they lie inside the domain of the rapid-change regime; there is no 
crossover line going through this point. This result is in agreement with the analysis of 
[T8] and in disagreement with [9] [H]. Possible explanation of the discrepancy between 
existing results was proposed in [19] : it was argued that existence of a crossover line going 
through the Kolmogorov point depends on the specific choice of the model parameters: 
if the amplitudes in the velocity correlators are related to the IR scale, the crossover 
disappears. The possibility that the stability domains of fixed points indeed change if 
the IR scale is introduced into the velocity correlators, was also discussed within the 
RG approach to the isotropic model; see section VIII in the e-print version of [23]. 

In this connection we stress that we found no crossover line, going through the 
Kolmogorov point, although in our approach the amplitudes in the velocity correlators 
(which play the part of the coupling constants) are related to the UV scale, see (12.191) and 
(19. 3p . In the terminology of [19], the correlators are generalized at the dissipation length. 
It is this choice that provides the agreement between the RG and other approaches to 
Kraichnan's model (for the discussion and comparison of various approaches, see |29j). 
In particular, if the UV scale A ~ l/l was replaced by the IR scale m ~ 1/C in 
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02.19p . the invariant charge f !7.5p in the inertial range would tend to zero instead of the 
nontrivial fixed point w* from (I4.8p . and the IR behaviour of the model would be trivial 
at all. Thus the problem requires further investigation. 

Two concluding remarks are in order. Although our results are exact, they are 
derived by RG and OPE resummations of the original perturbation series, and, in 
principle, their range of validity can be restricted by some boundaries in the £-77 plane, 
whose existence and location cannot be determined within the perturbative approach 
itself. The behaviour of the model can notably change, for example, for 77 > 2 (where 
the dispersion law becomes abnormal) or £ > 2 (where the eddy diffusivity becomes IR 
divergent). On the other hand, we mostly studied the structure functions, determined by 
the statistics of relative motion of the particles in the flow. The behaviour of individual 
particles is more subtle and sensitive to the details of the velocity statistics [T7] . 

In order to understand deeper the difference in the IR behaviour of the passive scalar 
in a weakly anisotropic velocity ensemble [20] and the shear flow of the type [9], it would 
be desirable to construct a more general model, which included them both as special 
limiting cases. Such a model is expected to demonstrate some kind of crossover between 
the IR behaviour described above for the latter case, and the anomalous (multi) scaling 
behaviour for the former one. It is not yet clear, however, how to do this. Probably, 
the idea [13] to approximate the isotropic case by a family of shear flows averaged with 
respect to their shearing directions will be useful here. This work remains for the future. 
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